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Abstract. We have analyzed the effects of photoelectric heating by dust grains in photoionization models of planetary neb- 
ulae. We have shown that this process is particularly important if planetary nebulae contain a population of small grains. The 
presence of such grains would solve a number of problems that have found no satisfactory solution so far: i) the thermal energy 
deficit in some objects inferred from tailored photoionization modelling; ii) the large negative temperature gradients inferred 
directly from spatially resolved observations and indirectly from integrated spectra in some planetary nebulae; iii) the fact that 
the temperatures derived from the Balmer jump are smaller than those derived from [O III] A4363/5007; iv) the fact that the 
observed intensities of [O I] A6300 are often larger than predicted by photoionization models. In presence of moderate density 
inhomogeneities, such as inferred from high resolution images of planetary nebulae, photoelectric heating would boost the 
temperature of the tenuous component, which would then better confine the clumps. The temperature structur e of s uch dusty 



and filamentary nebulae would solve the long standing problem of temperature fluctuations posed by Peimbert(1967) 
Ksy words. ISM: abundances; (ISM:) planetary nebulae: general ~ dust; Infrared: general; Radiative transfer. 



1. Introduction 

Infrared observations of planetary nebulae have established 
the presence of dust in these objects (see e.g. reviews by 
Barlow 1983 and Roche 1989). Infrared spectra of planetary 



nebulae clearly show continuum emission in the mid- to far- 
infrared range due to dust thermal emission (see e.g. Pottasch et 
al. 1984). Superimposed on this continuum are dust features at- 

Waters 



tributed to silicate dust particles (e.g. Aitken et al 
et al 



(e.g. Gillett et al.l973 



1979 



1998|, Cohe n et a l. 1999() or ca rbon-b ased dust pa rticles 

Cohen et al 



1998 



1999 



Waters et al, 

Szczerba et al. |2001| ). The large depletion factor s of refractory 
elements such as Fe, Mg and Si (Shi elds 1978 Pequignot & 



Stasiris ka|1980| Shields |1983[ Barlow |1993| ) and Ca (Aller & 
Czyzak 1983 Volk et al. 1997 ) indicate that dust is intimately 
mixed with the ionized gas. 

In a former paper (Stasiriska & Szczerba 1999), we per- 



formed a statistical analysis of a large sample of planetary neb- 
ulae with available IRAS photometric data by comparison with 
a grid of photoionization models for dusty planetary nebulae in 
which we computed the infra-red emission of dust heated by 
the absorbed stellar photons. The conclusion was that plane- 
tary nebulae contain dust in a proportion from about 10^^ up to 
10~^ by mas s (ado pting the canonical grain size distribution of 
Mathis et al.|l977|- hereafter the MRN distribution), and that 



there is no sign of an evolution of the dust content on a time 
scale of about Iff* yrs, contrary to earlier suggestions (Natta & 
Panagia|l98ll Pottaschll984l Lenzuni et al. |1989|). 
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That dust in ionized plasmas plays an important role in 
the thermal balance has been known since Spitzer(1948), 
who layed the basis for a quantitative modelling of the 
effects of heating and cooling of dusty plasmas. Baldwin 
et al. (|1991|), Borkowski & Harrington (|l 99 Binette et 



al. (1993), Viegas & Contini ( 1994i and more recently Dopita 
& Sutherland (2000) have included the effects of dust on the 



thermal balance of the gas in their photoionization codes. Of 
these, only the studies of Borkowski & HaiTington(1991) and 
Dopita & Sutherland ( 2000| ) deal explicitly with planetary neb- 
ulae. Borkowski & HaiTington(1991) studied a very specific 
object, the galactic halo planetary nebula IRAS 18333—2357 
located in the globular cluster M22, which is hydrogen-poor 
and extremely dusty. They showed that in this object, the heat- 
ing of the nebular gas is essentially provided by the photo- 
electric emission from dust grains. Dopita & Sutherland ( 200C| ) 
studied photoelectric heating by dust in normal planetary neb- 
ulae. They found its role to be important especially if one as- 
sumes the presence of small grains in a much higher proportion 
than inferred from the canonical MRN distribution. 

In this paper, we develop further on the ideas of Dopita 
& Sutherland (2000), addressing the temperature distribution 
inside planetary nebulae in more detail. This study is moti- 
vated by the fact that there are a number of unsolved prob- 
lems relating to the temperature structure of planetary nebulae. 
The best known concerns the temperature fluctuations, whose 



existence has been postulated by Peimbert(1967) to explain 



the difference between the temperature T^[oiu] derived from 
the [O III] A4363/5007 line ratio and the one derived from the 
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Balmer jump. Using Peimbert's formulation, one derives a tem- 
perature fluctuation typically of 0.03-0.04 in planetary neb- 
ulae (Liu & Danziger 1993[ , see also Esteban 2001 for a re- 
cent review). Although the Peimbert's approach relies on sim- 
plifications which are not necessarily justified (e.g. Kingdon 
& Ferland 1995) and the real temperature variations may not 
be adequately described by his scheme (Stasinska 2001 ), the 
mere fact that the Balmer jump temperature is significantly 
lower than Ti loiii] in ^ large number of planetary nebulae in- 
dicates that their temperature is far from uniform. While pho- 
toionization models do predict temperature gradients in plan- 
etary nebulae, as a result of the gradual modification of the 
ionizing radiation field and of the radial distributions of the 
ions responsible for the heating and the cooling of the gas, the 
computed gradient is too small to account for the difference 
between the Balmer jump temperature and T^[om] - Moderate 
density fluctuations are unable to produce temperature fluctua- 
tions on a t'^ = 0.03 level (Kingdon & Ferland 1995). Very high 



density clumps (n ^ 10 cm^-^) can bias upwards the tempera- 
ture derived from the [O III] A4363/5007 ratio due to collisional 
deexcitation of the [Olll] A5007 line (Viegas & Clegg |l994 - 
However, the existence of such high densities in the 0++ emit- 
ting zone is not supported by the density sensitive line ra- 
tios observed in planetary nebulae both from the infrared fine 
structure [Olll] lines and from the [Ariv] A4741/4713 dou- 



blet (e.g. Liu et al. 2001b). The magnitude of the temperature 
gradient observed in planetary nebulae either indirectly from 
a comparison of T^-[om] with Tj-jnh]' the temperature derived 
from [Nil] A 5755/6 584 (see e.g. Kingsburgh & Barlow |1 994 
or Pena et al. 2001) or directly from spatially resolved obser- 
vations (Liu et al. 200C, Garnett & Dinerstein 2001) is often 



difficult to account for by photoionization models. 

In this paper, we explore whether photoionization models 
including the effects of dust on the thermal balance of the emit- 
ting gas can answer some of these questions. 

In the next section, we briefly describe the photoionization 
code and dust properties. In Section 3, we study the tempera- 
ture profiles of planetary nebulae with smooth density distri- 
butions. In Section 4, we discuss the temperature profiles of 
planetary nebulae with small scale density inhomogeneities. In 
Section 5, we discuss the effects of dust on methods of analysis 
of planetary nebulae based on emission line ratios. In Section 6 
we discuss the temperature fluctuation parameter. Concluding 
remarks are given in Section 7. 

2. The building of dusty pliotoionization models 

2. 1. The treatment of dust in the photoionization 
code 

The present paper is based on results from the photoioniza- 
tion code PHOT O using the same atomic data as Stasinska & 
Leitherer(1996). All the computations have been done under 
the assumption of spherical symmetry and adopting the out- 
ward only approximation for the diffuse ionizing radiation. 
In the version of PHOTO used by Stasinska & 



grains, and their conversion into thermal radiation from the 
grains. For the present paper, we have also incorporated the 
effects of grains on the thermal balance of the gas, by imple- 
menting the physical processes between dust, gas and radiation 
relevant for dusty ionized regions as described in Appendix C 
of Baldwin et al. (1991) - except that we do not consider an 



average grain size like former authors, but consider all the par- 
ticle sizes explicitly. In short, the steady state grain potential 
is determined self-consistently for each dust particle size and 
the corresponding heating of the electron plasma by the photo- 
electric emission from dust is computed. As a check, the same 
treatment of dust has been applied to the photoionization code 



written by Szczerba (see e.g. Gfsicki et al. 199(:)- Both codes 



give consistent results and reproduce the Baldwin et al. ( |1991| ) 
model for the Orion Nebula satisfactorily (their Fig. 16). 

For simplicity, in this p aper we consider only graphite 
grains (Draine & Laor 1993 ) in our planetary nebulae mod- 
els. We include two population of grains. One is composed of 
graphite grains with a bulk density of 2.5 g cm^'' and sizes dis- 
tributed between 0.01 and 0.25 /im according to the classical 
MRN model (hereafter referred to as standard grains). The sec- 
ond one consists of small grains, as introduced by Dopita & 
Sutherland (200C). We assume that the small grain have a bulk 
density of 1 g cm^'^ and that their sizes are distributed accord- 
ing to the same MRN law but between 0.001 and 0.01 /im. The 
formulae describing the photoelectron energy distribution and 



the yield are from Dopita & Sutherland ( 200C ). As explained by 



these authors, in comparison with the formulation of Baldwin 



et al. (1991) the adopted approach seems to be more realistic. 



Szczerba (1999), we had adopted the formulation of Harrington 



et al. (1988) to compute the absorption of photons by dust 



2.2. The planetary nebula models 

We consider two families of photoionization models for plan- 
etary nebulae. The first one corresponds to a rather dense 
nebula ionized by a star that is not very hot. The reference 
model has a density n(H) =10** cm^^^ and a stellar temperature 
T*=5 10'' K. The second one corresponds to a more diluted 
nebula ionized by a hotter star We have chosen for the refer- 
ence model n(H) =5 lO'^ cm^'^ and = 10^ K. The stars are 
assumed to radiate as blackbodies. The inner radius of the neb- 
ula is chosen to be 5 10^® cm for the former case, and 8 10^® cm 
for the latter one. In both cases, we assume that the total num- 
ber of H Lyman continuum photons emitted by the star, Q{}i^) 
is equal to 3 10**^ phs^^ and that the nebulae are ionization 
bounded (the calculations are stopped when the proportion of 
neutral hydrogen exceeds 0.03). This implies that the total neb- 
ular mass is larger than ~ 0.1 AIq for the first case, and larger 
than ~ 0.3 Mq for the second case, which are reasonable val- 
ues. Actually, the local properties of our models do not depend 
on whether the nebula is ionization- or density-bounded, but 
the integrated properties do depend on that, since the low ex- 
citation parts are suppressed for density-bounded nebulae. In 
Table 1 and Table 2 we have collected the key parameters used 
in the models computed by us for =5 10** K and 10^ K, re- 
spectively. These tables allow to easily relate the name of the 
model to the physical parameters adopted in that model. 
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Table 1. Definition of the models with = 5 lO'^ K. 



model 


Aa 


Ba 


Ca 


Da 


Ea 


Fa 


Ga 


T* [K] 


510* 


510* 


5 10* 


5 10* 


5 10* 


5 10* 


510* 


Q(H°) [phs^M 


3 10*" 


3 10*" 


310*" 


310*" 


3 10*" 


310*" 


3 10*" 


-Rin [cm] 


5 lO^'^ 


5 lO^^ 


510" 


510" 


5 10" 


510" 


510" 


n(H) [cm-3] 


10* 


10* 


10* 


10* 


Eq (3) 


Eq (4) 


Eq (4) 


dust properties 


no grains 


standard grains 


standard + 


standard + 


no grains 


no grains 


standard + 








small grains 


small grains 






small grains 


md/mn 


0. 


10"^ 


5 10^3+5 10-3 


5 10-*+5 10-* 


0. 


0. 


5 10"3+5 10-3 


Table 2. Definition of the models with T, = 10 












model 


Ab 


Bb 


Cb 


Db 


Eb 


Fb 


Gb 


T. [K] 


lO'^ 


10=' 


10' 




lO'' 


10=' 


10^ 


0(H°)[phs-i] 


3 10*" 


3 10*" 


3 10*" 


3 10*" 


3 10*" 


3 10*" 


310*" 


-Rin [cm] 


810^^ 


810" 


810" 


810" 


810" 


810" 


8 10" 


7i(H) [cm-3] 


510=* 


510^ 


510^ 


5103 


Eq (3) 


Eq (4) 


Eq (4) 


dust properties 


no grains 


standard grains 


standard + 


standard + 


no grains 


no grains 


standard + 








small grains 


small grains 






small grains 


md/rriH 


0. 




5 10-3+5 10^3 


5 10-*+5 10"* 


0. 


0. 


5 10^3+5 10-3 



Our models have the following chemical composition: 
He/H = 0.115; C/H = 5.50 10"*; N/H = 2.24 10"*; O/H 
= 4.7910-*; Ne/H = 1.2310"*; Mg/H = 3.80 10"^; Si/H 
= 3.5510-5; s/H = 8.3210-6; ci/H = 1.86 IQ-"; Ar/H = 
2.46 10-6; Fe/H = 4.68 10-'', which corresponds to the average 
chemical composition of th e Gal actic planetary nebulae from 
the Kingsburgh & Barlow ( |1994 sample (Mg, Si, CI and Fe 
have solar abundances). This chemical composition is assumed 
throughout or study, regardless of the dust abundance. The rea- 
son for adopting such a policy is its simplicity. Of course, 
if metallic or carbon atoms are condensed into grains, their 
gas phase abundances are correspondingly depleted. However, 
such elements as Si, Mg, Fe do not play a dominant role in the 
energy balance. Carbon, on the other hand, can play an impor- 
tant role, especially in carbon-rich planetary nebulae. However, 
we are not interested here to investigate what would happen to 
a planetary nebula of given chemical composition if part of its 
carbon content were locked into grains. We rather wish to ex- 
plore what are the effects of dust on the temperature structure 
of planetary nebulae. 

In Tables || and ^ we summarize some interesting features 
of the models discussed in this paper We list the total H/3 
fluxes in ergcm--^ s-^ and the fluxes in the four IRAS bands 
at 12, 25, 60 and 100 /im, assuming that the nebulae lie at a 
distance of 1 kpc. Below we give the intensities relative to H/3 
of heUum and oxygen lines. These are sufficient to give a fair 
understanding of the behaviour of the emission line spectra. 
We also list the intensities of the most important lines of C, 
N and Ne used in abundance determinations. We do not com- 
pute the reddening by internal dust, so that the listed line in- 
tensities are relevant to observations corrected for reddening 



(the dereddening procedure, ideally consisting in fitting the ob- 
served Balmer decrement to the theoretical one, actually pro- 
duces a small error in the ratios of lines arising from differ- 
ent ions). Then we list the two temperature sensitive fine ra- 
tios [O III] A4363/5007 and [N ll] A5755/6584, and the Balmer 
jump defined as (FA3642-FA3648)/F(H/3), where FA3642 
and f A3648 are the nebular continuum fluxes at 3642 and 
3648 A. Below are given the two density sensitive line ratios 
[S II] A6731/6717 and [Ariv] A4741/4713. 

Next we give the values of the average ionic temperatures 
ro(H+), ro(0+) and To(0++), where ro(X+^) is defined as: 



To(X+-') = 



J n(X+J)nc dV 



(1) 



and To is the local electron temperature, tIq the electron density, 
n(X+^ ) stands for the ion concentration, the integration being 
performed over the total volume. 

The last three lines give the values of t^(H+), t^(0+) and 
t'^iO++), where t'^{X+^) is defined as: 



t2(X+J') 



To(X+J')2 /n(X+J>edl/ ■ 



(2) 



3. The effect of dust on the temperature profiles 
of planetary nebulae with smooth density 
distributions 

We first investigate the simplest case of a uniform density 
distribution. Figure ^ shows the results of dust-free models 
(Models Aa and Ab) for the n(H) = 10* cm-^, =5 10* K 
case (upper panels) and for the n(H) =510'^ cm-^, T* = 10^ K 
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Table 3. Selected properties for models with = 5 10'' K. 



model Aa model Ba model Ca model Da model Ea model Fa model Ga 





1.149E-09 


8.082E-10 


1.774E-10 


7.250E-10 


1.137E-09 


1.145E-09 


2.202E-10 


TpIR.AS 
^12 
rilRAS 
^25 
rfRAS 
^60 
piIRAS 


6.175E+00 
5.667E+00 
1.473E+01 
3.302E+00 


2.819E+01 
3.885E+02 
3.865E+02 
7.361E+01 


1.420E+02 
6.182E+02 
3.445E+02 
5.499E+01 


5.914E+01 
2.484E+02 
1.376E+02 
2.173E+01 


5.155E+00 
4.572E+00 
1.048E+01 
2.200E+00 


5.342E+00 
4.603E+00 
1.007E+01 
2.290E+00 


1.400E+02 
5.926E+02 
3.154E+02 
4.912E+01 


He I A5 876 
He II A4686 


1.680E-01 
1.728E-03 


1.711E-01 
2.405E-03 


1.890E-01 
8.745E-03 


1.745E-01 
2.556E-03 


1.693E-01 
1.763E-03 


1.689E-01 
1.730E-03 


1.882E-01 
7.367E-03 


[Oi] A6300 
[On] A3726,3729 
[On] A7320,7330 
On A4651 
III] A 1663 
[Oiii] A4363 
[Oiii] A5007 
[O III] A52 fim 
[Oiii] A88/im 


6.771E-03 
5.428E-01 
6.371E-02 
4.101E-03 
9.149E-03 
1.530E-02 
4.627E+00 
3.486E-01 
4.392E-02 


8.911E-03 
6.018E-01 
7.428E-02 
4.165E-03 
1.385E-02 
2.067E-02 
5.410E+00 
3.650E-01 
4.603E-02 


6.110E-02 
1.466E+00 
2.564E-01 
4.302E-03 
6.448E-01 
2.942E-01 
1.882E+01 
5.100E-01 
6.481E-02 


1.615E-02 
7.069E-01 
9.629E-02 
4.396E-03 
3.619E-02 
4.051E-02 
7.625E+00 
4.107E-01 
5.189E-02 


2.651E-03 
4.573E-01 
6.211E-02 
4.450E-03 
1.134E-02 
1.870E-02 
5.217E+00 
2.502E-01 
2.936E-02 


1.350E-02 
4.612E-01 
1.053E-01 
3.649E-03 
9.393E-03 
1.553E-02 
4.264E+00 
2.380E-01 
3.058E-02 


7.165E-02 
1.187E+00 
3.362E-01 
3.779E-03 
4.560E-01 
2.148E-01 
1.461E+01 
3.276E-01 
4.259E-02 


Cm] A 1909 


3.095E-01 


4.357E-01 


9.567E+00 


9.298E-01 


3.663E-01 


3.244E-01 


6.858E+00 


[Nil] A6584 


1.334E+00 


1.391E+00 


2.514E+00 


1.466E+00 


1.071E+00 


1.589E+00 


2.790E+00 


[Neiii]A3869 


3.679E-01 


4.486E-01 


2.119E+00 


6.829E-01 


4.192E-01 


3.584E-01 


1.689E+00 


[Oiii] A4363/5007 
[Nil] A5755/6584 
(i^A3642-FA3648)/F(H/3) 
[S II] A6731/6717 
[ArlV] A4741/4713 


3.308E-03 
1.132E-02 
4.902E-03 
1.991E+00 
1.691E+00 


3.822E-03 
1.215E-02 
4.801E-03 
1.980E+00 
1.673E+00 


1.564E-02 
1.940E-02 
4.026E-03 
1.860E+00 
1.500E+00 


5.313E-03 
1.407E-02 
4.593E-03 
1.924E+00 
1.629E+00 


3.585E-03 
1.248E-02 
4.852E-03 
1.896E+00 
2.340E+00 


3.643E-03 
1.752E-02 
4.855E-03 
2.097E+00 
1.885E+00 


1.470E-02 
2.535E-02 
4.132E-03 
1.986E+00 
1.598E+00 


ro(H+) 
ro(0+) 

ro(o++) 


8151. 
7972. 
8232. 


8427 
8184. 
8529. 


11807. 
10010. 
12360. 


9070. 
8672. 
9204. 


8279. 
8116. 
8334. 


8268. 
8162. 
8335. 


11208. 
9870. 
11836. 


t^(H+) 
t^(0+) 
t'^{0++) 


8.0E-4 
1.3E-3 
2.9E-4 


1.3E-3 
1.6E-3 
6.7E-4 


4.6E-2 
1.2E-2 
4.3E-2 


4.5E-3 
1.9E-3 
4.4E-3 


7.0E-4 
1.6E-4 
2.4E-4 


6.6E-4 
7.7E-4 
4.3E-4 


4.4E-2 
7.3E-3 
4.6E-2 



case (lower panels). The left panels show the electron temper- 
ature (thick line) and hydrogen number density (thin line) as a 
function of the distance R to the star. The right panels show the 
ionization fraction of He+ (thick line) and 0++ (thin line). 

In Figure^, we display with the same conventions the re- 
sults of dusty models (Ba and Bb) in which dust is composed 
of graphite grains with the MRN size distribution between 0.01 
and 0.25 /^m and a dust-to-hydrogen mass ratio, rnd/mn, of 
10~^. Clearly, the size of the ionized zone is shrunken with 
respect to the dust-free models shown in Fig.[l|, because with 
such md/mn the grains significantly compete with gas for the 
absorption of the Lyman continuum photons. However, there 
is no strong effect neither on the temperature structure of the 
nebula nor on the ratios of emission lines with respect to H/3, 
as can be seen in Table |] and ^. The only notable difference 
is in the H/3 and the IRAS fluxes. The H/3 flux is reduced in 



the dusty models by about 30%, while the fluxes in the IRAS 
bands are increased by large amounts due to the thermal emis- 
sion of the dust grains that have been heated by the stellar and 
nebular photons (note that in the dust-free models, the fluxes in 
the IRAS bands are not zero due to the contribution of the fine 
structure lines from atomic ions). 

Figure|3] shows models Ca and Cb which differ from mod- 
els Ba and Bb only in that half of the dust mass is now in small 
grains of size from 0.001 to 0.01 /im. That is to say, large grains 
and small grains coexist, with a dust-to-hydrogen mass ratio 
of 5 10^'^ for each component. Now, the H/3 fluxes are con- 
siderably reduced, by a factor 4 in the = 10^ K model and 
by a factor 6 in the = 5 lO"* K one, with respect to the dust- 
free models. Indeed, the mass absorption efficiency of small 
graphite grains is much larger than that of large ones (see 



Draine 8c Laorl993). In addition, small grains are heated to 
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Table 4. Selected properties of models with = 10^ K. 



5 



model Ab model Bb model Cb model Db model Eb model Fb model Gb 





1.185E-09 


8.961E-10 


3.204E-10 


8.930E-10 


1.163E-09 


1.169E-09 


2.939E-10 


TpIR.AS 
^12 
rilRAS 
^25 
rfRAS 
^60 
piIRAS 


1.064E+01 
1.879E+01 
3.721E+01 
1.061E+01 


1.162E+01 
2.224E+02 
4.668E+02 
1.199E+02 


4.887E+01 
5.616E+02 
5.697E+02 
1.185E+02 


2.260E+01 
1.917E+02 
1.936E+02 
4.045E+01 


1.057E+01 
1.193E+01 
3.154E+01 
8.651E+00 


8.337E+00 
1.524E+01 
2.452E+01 
6.967E+00 


8.595E+01 
5.785E+02 
4.604E+02 
8.842E+01 


He I A5 876 
He II A4686 


1.636E-01 
1.082E-01 


1.606E-01 
1.362E-01 


1.432E-01 
3.212E-01 


1.616E-01 
1.385E-01 


1.635E-01 
1.123E-01 


1.672E-01 
1.109E-01 


1.471E-01 
3.363E-01 


[Oi] A6300 
[On] A3726,3729 
[On] A7320,7330 
On A4651 
III] A 1663 
[Oiii] A4363 
[Oiii] A5007 
[O III] A52 fim 
[Oiii] A88/im 


5.159E-02 
l.UOE+OO 
1.086E-01 
4.830E-03 
1.059E-01 
9.251E-02 
1.255E+01 
8.893E-01 
1.379E-01 


5.559E-02 
1.189E+00 
1.153E-01 
4.670E-03 
1.206E-01 
1.005E-01 
1.284E+01 
8.716E-01 
1.353E-01 


1.044E-01 
1.982E+00 
2.083E-01 
3.905E-03 
4.862E-01 
2.578E-01 
1.859E+01 
8.428E-01 
1.321E-01 


6.262E-02 
1.306E+00 
1.315E-01 
4.677E-03 
1.638E-01 
1.262E-01 
1.445E+01 
9.008E-01 
1.403E-01 


2.313E-02 
8.950E-01 
7.963E-02 
5.192E-03 
1.365E-01 
1.144E-01 
1.439E+01 
7.739E-01 
1.143E-01 


7.552E-02 
1.022E+00 
1.821E-01 
4.557E-03 
1.043E-01 
9.109E-02 
1.189E+01 
5.950E-01 
9.172E-02 


1.615E-01 
1.757E+00 
3.282E-01 
3.638E-03 
4.718E-01 
2.475E-01 
1.727E+01 
5.527E-01 
8.593E-02 


Cm] A 1909 


2.093E+00 


2.289E+00 


6.975E+00 


2.931E+00 


2.579E+00 


2.095E+00 


6.675E+00 


[Nil] A6584 


1.829E+00 


1.955E+00 


3.042E+00 


2.086E+00 


1.217E+00 


2.272E+00 


3.632E+00 


[Neiii]A3869 


1.341E+00 


1.403E+00 


2.355E+00 


1.609E+00 


1.482E+00 


1.353E+00 


2.307E+00 


[Oiii] A4363/5007 
[Nil] A5755/6584 
(i^A3642-FA3648)/F(H/3) 
[S II] A6731/6717 
[ArlV] A4741/4713 


7.374E-03 
1.688E-02 
4.279E-03 
1.690E+00 
1.183E+00 


7.828E-03 
1.668E-02 
4.256E-03 
1.689E+00 
1.179E+00 


1.387E-02 
1.877E-02 
3.965E-03 
1.674E+00 
1.147E+00 


8.729E-03 
1.764E-02 
4.165E-03 
1.685E+00 
1.172E+00 


7.948E-03 
1.826E-02 
4.228E-03 
1.446E+00 
1.430E+00 


7.660E-03 
2.123E-02 
4.269E-03 
1.809E+00 
1.435E+00 


1.433E-02 
2.311E-02 
3.979E-03 
1.828E+00 
1.406E+00 


ro(H+) 
ro(0+) 

ro(o++) 


10410. 
10097. 
10310. 


10602. 
10045. 
10489. 


12906. 
10575 
12576. 


10988. 
10312. 
10874. 


10615. 
10566. 
10523. 


10462. 
10164. 
10347. 


13120. 
10661. 
12597. 


t^(H+) 
t^(0+) 
t'^{0++) 


4.9E-3 
2.4E-3 
2.1E-3 


6.7E-3 
2.9E-3 
2.6E-3 


3.4E-2 
8.1E-3 
1.4E-2 


7.3E-3 
3.2E-3 
2.7E-3 


3.7E-3 
2.7E-3 
2.4E-3 


5.9E-3 
1.3E-3 
2.6E-3 


6.5E-2 
7.0E-3 
1.6E-2 



larger temperatures, because the heating of grains is propor- 
tional to their surface while the cooling is proportional to their 
volume. This translates into much larger fluxes in the 12 /xm 
band (by a factor of about 4) and in the 25 fim band (by a fac- 
tor of about 2) when comparing with models computed solely 
with the MRN dust, as seen in Tables]^ and ^ 

What is striking from the comparison of Fig.jl with Figs. |l] 
and ^ is the very significant increase in the gas temperature 
when a population of small grains is included. Figure^ shows 
the contribution of the gas heating due to the photoelectric ef- 
fect with respect to the total heating (thick line) and the ratio 
between gas cooling due to grain-gas collisions and the total 
cooling (thin line) as a function of the distance to the cen- 
tral star As expected, the contribution of photoelectric heat- 
ing is largest close to the star Indeed, the heating due to 
photoelectrons from dust grains per unit volume of gas, Fd, 



is proportional to the number density of dust grains and to 
the intensity of the stellar radiation field. The heating due to 
photoionization of hydrogen, Fh, is proportional to the num- 
ber density of neutral hydrogen and to the stellar radiation 
field. Using the ionization equilibrium equation for hydrogen, 
one then finds that Fh is simply proportional to the square 
of the gas density. It follows that the ratio Fd/Fn is sim- 
ply proportional to {nid/m,^) U, where U is the local ion- 
ization parameter, defined as (5(H°)/(4 tt i?^ rio c) (c is the 
speed of light). Therefore, the photoelectric effect becomes in- 
significant in fully ionized plasma at large distances from the 
central stars. This results in an important electron tempera- 
ture gradient, which was also seen in the models of Dopita 
& Sutherland ( ^000[ ). Consequently, [O III] A4363/5007 is in- 
creased, and the intensities of coUisionally excited lines such as 
[O III] A5007, [O III] A4363 or O III] A 1663 are increased with 
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Fig. 1. The temperature and ionization structure of dust-free and constant density model planetary nebulae (Models Aa and Ab). 
Upper panels: models with n(H) = 10'* cm^'' and = 5 10'' K. Lower panels, models with ri(H) =5 10"^ cm~^ and T» = 10^ K. 
Left: electron temperature (thick line) and hydrogen concentration (thin line) as a function of the distance to the star; right: 
ionization fraction of He+ (thick line) and 0++ (thin line). Note that here as well as in the other figures of the paper, the 
horizontal scale matches the size of the ionized region and is therefore different for each model, while the vertical scale is always 
the same. 



respect to H/3, as seen in Tabled and Fi As already noted by 
previous workers (see Baldwin et al. |l991 ) cooling due to gas- 
grain collisions can also be important in the energy budget of 
the gas, and the net effect of the grains on the gas temperature 
is not necessarily positive. In the models shown in Fig. ^ the 
net effect is slightly negative in the zone. 

Another noteworthy effect of the presence of small grains is 
the strong heating of the gas due to the photoelectric effect on 
grains near the ionization front. In this region, the Lyman con- 
tinuum photons have been exhausted, and the only photons that 
are present are stellar and nebular photons below the Lyman 
limit (including the Ly a radiation) which are not absorbed by 
hydrogen, but can heat the gas via the photoelectric effect from 
dust grains. This results in an enhancement of the [O l] A6300 
line, as seen in Tables 3 and 4. 

The heating efficiency of the small dust grains is such that 
even if present in a very small amount, they can affect the gas 
thermal balance significantly. We have run models Da and Db 
which are identical to models Ca and Cb except that the dust-to- 
hydrogen mass ratio of each grain population is 5 lO^"' instead 
of 5 10-3. Table I shows that, in the = 5 lO"* K models, such 
a small amont of dust increases the [O III] A4363/5007 ratio 
from 3.31 lO^^^ (dust free case) to 5.31 lO^^^ and the electron 
temperature ro(0++) from 8200 to 9200K. Such a difference 
is quite significant if one is interested in an accurate model fit- 
ting of a nebula. When comparing the models with T^, = lO'^ K 



in Table 4 (models Ab and Db), the results are less spectacu- 
lar: the increase in To(0++) is only by 500 K. This is because 
the photoelectric heating contribution is important in the in- 
ner regions, where oxygen is in the form of 0+^+ and not of 
The effect of dust heating is better seen on lines of more 
charged ions, such as C IVA1550 whose intensity with respect 
to H/3 (not shown in Table ^ is increased by a factor 2. 

One way to measure a temperature gradient in a nebula is 
to compare ro(0++) and To(0+). We see from Table 3 and 4 
that models with small grains increase the To(0++)/To(0+) 
ratio. As a matter of fact, as already stated by Dopita & 
Sutherland (2000), dust could help explaining the observed 



temperature gradients in planetary nebulae. Indeed, the ratio 
ri.[oiii]/T'r[Nn] takes values between 1 and 1.5 in the planetary 
nebulae sample of Kingsburgh & Barlow (1994). Such gradi- 
ents are very difficult to produce with dust-free photoionization 
models. The most natural cause to investigate is density gradi- 
ents, since the temperature is expected to be higher in regions 
of higher density, due to collisional deexcitation of the cooling 
lines. Figure 5 displays such models (Ea and Eb), in which the 
density at given distance R from the star is given by: 



n(H) = ?T.(H)i,i exp ■ 



R — i?i] 



(3) 



where is the inner radius of the nebula and h the thickness 
scale of the envelope chosen to be 5 10*^ cm for T, = 5 10'' K 
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Fig. 2. The temperature and ionization structure of the same constant density models as in Fig. |l] but containing graphite grains 
with the MRN size distribution between 0.01 and 0.25 /im and a dust-to-hydrogen mass ratio of 10^^ (Models Ba and Bb). The 
presentation is the same as for the models shown in Fig. Ill 




Fig. 3. The temperature and ionization structure of the same constant density models as in Fig. ^ but in which there are two 
populations of grains: small grains and large grains, each population having a dust-to-hydrogen mass ratio of 5 10^'^ (Models Ca 
and Cb). The presentation is the same as for the models shown in Fig. ^ 
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Fig. 5. The temperature and ionization structure of dust-free models with a density gradient, models Ea and Eb (see text). The 
presentation is the same as for the models shown in Fig. Ill 



and 10^^ cm for = 10^ K, respectively. These models have 
roughly the same outer radius and H/3 luminosity as mod- 
els Aa and Ab. Our computed models show that the effect 
of such a density gradient is far from sufficient to produce 
an appreciable temperature gradient. Of course, one can think 
of much steeper density gradients, but in order to model a 
given object, the gradient should be compatible with the ob- 
served surface brightness distribution and the density derived 
from [Sll] A6731/6717 and [Ariv] A4741/4713. Our experi- 
ence is that it is generally very difficult to reproduce the ob- 
served temperature gradient with a density gradient compatible 
with the observations. Another possibility is to invoke deple- 
tion or chemical abundance gradients. But then it is the outer 
parts of the nebulae which should be richer in heavy elements, 
which makes such an assumption rather unlikely. The hypothe- 
sis of photoelectric heating by small dust grains provides a nat- 
ural explanation to the outwards decrease of electron tempera- 
tures that is inferred from the observations. Of course, given all 
the unknowns as to the nature and physical properties of such 
grains, a detailed model fitting including small dust grains ob- 
viously would bear important uncertainties. 

4. The effect of dust on the temperature profiles 
of planetary nebulae with filamentary 
structures 

As mentioned in the introduction, density condensations have 
been proposed to explain the temperature fluctuations inferred 
from the differences between the temperatures derived from 
various indicators. But in order to significantly enhance the 



temperature derived from [O III] A4363/5007, very high den- 
sities (~ 10® cm~^) are required. Such densities are excluded 
by observations (e.g. Liu et al. 2001a). Modest density fluc- 
tuations, though, are obviously present, as seen on numerous 
high resolution images of planetary nebulae which reveal knots 
and filaments. In a dust-free nebula, the filaments will have a 
lowered ionization compared to the ambient gas, but the tem- 
perature will not be largely different from that of the ambient 
medium. We have modeled such a situation by assuming that 
the filaments are concentric shells and that the density distribu- 
tion in the nebula is given by: 

n(H) = n(H),„in + <5n(H) (cos^^ (4) 

where q is taken equal to 15 in order to produce filaments that 
occupy a small volume of the nebula. 

Figure^ shows such filamentary model nebulae (models Fa 
and Fb) whose density distribution has been chosen to provide 
the same H/3 luminosity and size as models Aa and Ab, with a 
density contrast (5n(H) of a factor 5. Note that such completely 
ionized filaments are far from pressure equilibrium with the 
ambient gas, so their lifetime is very short (about 30 years for 
filaments of thickness lO^^cm) which means that they contin- 
uously form and disappear during the lifetime of the planetary 
nebula. If however the gas contains a population of small dust 
grains, the temperature inside the filaments is expected to be 
smaller than that of the ambient gas, since the importance of 
gas heating due to the photoelectric effect is roughly propor- 
tional to the ionization parameter, as stated above. Therefore, a 
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R/IO^^ [cm] R/10^^ [cm] 

Fig. 6. The temperature and ionization structure of the dust-free filamentary models Fa and Fb, which have the same total H/3 
flux as models Aa and Ab (see text). Same presentation as in Fig. |l| 




[cm] R/^0^' [cm] 

Fig. 7. The temperature and ionization structure of filamentary models (Ga and Gb) with similar input parameters as models Fa 
and Fb except that they contain dust with the same characteristics as models Ca and Cb. Same presentation as in Fig.[l|. 

decrease of U in the denser regions will result in smaller heat- the same input values as models Fa and Fb except that they 
ing due to the photoelectrons from dust grains. Figure^ illus- contain dust with the same characteristics as models Ca and 
trates such a case, by showing models Ga and Gb, which have Cb, respectively (i.e. small and large grains, each population 
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Fig. 4. The contribution of the heating due to the photoelectric 
effect from dust grains with respect to the total heating (thick 
line) and the proportion of cooling by grain-gas collisions with 
respect to the total cooling (thin line) as a function of the dis- 
tance to the central star for the models shown in Fig. I (Models 
Ca - top and Cb - bottom). 

with a dust-to-hydrogen mass ratio of 5 lO^''). Of course, the 
H/3 luminosity and the size of the ionized region is smaller than 
in models Fa and Fb, due to the fact that dust grains compete 
with hydrogen to absorb the H Lyman continuum photons. We 
clearly see the predicted effect: the temperature inside the fila- 
ments is lower than in the ambient medium. The magnitude of 
the difference decreases with the distance to the ionizing star. 
Note that in such a model, the pressure discontinuities at the 
borders of the filaments are much smaller than in a dust-free 
model. In other words, filaments in a nebula containing small 
grains can survive longer than in a dust-free nebula. 

5. The effect of dust on the emission line 
analysis of planetary nebulae 

It is instructive to analyze our models with techniques used for 
the interpretation of the spectra of real planetary nebulae. The 
virtue of such "numerical observations" is that they eliminate 
the uncertainties in the atomic data, since they can be done us- 
ing exactly the same atomic data as used in the computations of 
the models. Such an approach has been applied in different con- 



texts by Gruenwald & Viegas ( 1998 ) and Perinotto et al. ( 1998 ) 



5. 1. Plasma diagnostics and abundance 
determinations from integrated spectra 

We have used the usual procedure of deriving Ti.[oiii] from 
[O III] A4363/5007, r,[Nii] from [Nil] A5755/6584, 7ie[sii] 
from [S II] A6731/6717 and n^[A,iv] from [Ariv] A4741/4713 
ratio. The emissivities of hydrogen and helium ions as well as 
those of doubly ionized species of C, N, O and Ne were com- 
puted using Ti.[oiii] and no[Ariv] while those of singly ion- 
ized species of C, N, O and Ne were computed using T,.[nii] 



and 



[SII]' 



The atomic abundance ratios were then derived 



from the ionic abundance ratios using the ionization correc- 



tion factors of Kingsburgh & Barlow ( |1994[ ). The abundance of 
O was obtained using the value of 0++/H+ as derived from 
[O III] A5007/H/?. The C/O, N/0 and Ne/0 ratios were derived 
from the Cm] A1909/[Oiii] A5007, [Nil] A6584/[Oii] A3727 
and [Nelll] A3869/[Olll] A5007, respectively. The results are 
displayed in Tables ^ and ||, respectively. For oxygen, we 
also show the determinations of 0++/H+ based on the other 
lines produced by 0"*"+ (including the recombination line 
O III A4651) but still using T^[om])- Since these lines have dif- 
ferent dependences on the electron temperature (and density), 
differences between these various determinations of the 0+"*" 
abundance in the same model from integrated "numerical spec- 
tra" are a good measure of the temperature variations in the 
models (a somewhat similar approach was adopted by Mathis 



Here we use our code ABELION to analyze the models com- 
puted with PHOTO with the same atomic data. 



et al. 1998) in the discussion of observational data on planetary 
nebulae. 

Inspection of Tables H and || reveals that the only cases 
where the computed O/H is significantly different from the 
input value are models Ca, Fa and Ga [} In the constant 
density case with small grains (model Ca), O/H is under- 
estimated by about 40%, because of the strong temperature 
gradient due to photoelectric heating by dust. In the dust- 
free case with filaments (case Fa), O/H is underestimated by 
30%, mainly because the determination of the 0+ abundance 
is biased by the presence of clumps of density higher than 
the critical density for the [O ll] A3727 line. The case where 
the difference between the derived O/H and the true one is 
the largest (almost a factor 2) is in model Ga, where small 
grains are present in a clumpy medium. Then, the temper- 
ature inhomogeneities are quite large and produce a signif- 
icant bias in the oxygen abundance determination. That the 
effect is mainly due to temperature inhomogeneities can be 
seen from the fact that O++(A1663)/O++(A5007) is 1.21 and 
O++(A4651)/O++(A5007) is 1.40, whereas these ratios are 
0.985 and 1.05, respectively, for the filamentary dust-free case. 
Note that in this model, the Balmer jump temperature (ob- 
tained by calibrating our theoretical Balmer jump index on 
isothermal photoionization models) is 10800 K, significantly 
smaller than the temperature derived from [O III] A4363/5007 
and [N ll] A5755/6584. As regards the models with = 10^ K, 
although qualitatively the effect of dust heating is similar to the 
one occuring in the models with =5 10"* K (see Figs.|| and 

' One does not expect the O/H determinations to be exactly equal to 
the input values, since the temperature characteristic of the emission 
of the H/3 line is not equal to Tr[o iii] and the correction for "unseen" 
oxygen ions introduces a slight error. 
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model Aa model Ba model Ca model Da model Ea model Fa model Ga 



-^r [O III] 


8242. 


8552. 


13504. 


9356. 


8330. 


8422. 


13145. 


^r[NII] 


8197. 


8434. 


10480. 


9079. 


8776. 


9051. 


11334. 


TlnfQ TTl 


9120. 


8837. 


6832. 


7444. 


6548. 


14782. 


10108. 


Tiof A r TA/"! 


11120. 


11091. 


10949. 


11012. 


20401. 


13830. 


12338. 


He+(A5876)/H+ 


l.UE-Ol 


1.13E-01 


1.07E-01 


1.14E-01 


1.12E-01 


1.12E-01 


1.08E-01 


He++(A4686)/H+ 


1.34E-04 


1.88E-04 


7.36E-04 


2.03E-04 


1.37E-04 


1.35E-04 


6.18E-04 


He/H 


1.12E-01 


1.14E-01 


1.08E-01 


1.15E-01 


1.12E-01 


1.12E-01 


1.08E-01 


0+(A3727)/H+ 


1.20E-04 


1.12E-04 


7.04E-05 


8.23E-05 


6.34E-05 


8.83E-05 


5.55E-05 


O++(A5007)/H+ 


3.26E-04 


3.30E-04 


2.61E-04 


3.32E-04 


3.57E-04 


2.77E-04 


2.18E-04 


0++(A1663)/H+ 


3.37E-04 


3.42E-04 


3.08E-04 


3.56E-04 


3.71E-04 


2.73E-04 


2.63E-04 


0++(A52 Atm)/H+ 


3.27E-04 


3.31E-04 


3.08E-04 


3.43E-04 


4.04E-04 


2.66E-04 


2.25E-04 


0++(A88 Atm)/H+ 


3.26E-04 


3.30E-04 


3.05E-04 


3.42E-04 


4.17E-04 


2.83E-04 


2.35E-04 


0++(A4651)/H+ 


3.28E-04 


3.33E-04 


3.48E-04 


3.53E-04 


3.56E-04 


2.92E-04 


3.06E-04 


0/H 


4.47E-04 


4.43E-04 


3.32E-04 


4.15E-04 


4.21E-04 


3.66E-04 


2.75E-04 



C/O 1.29E+00 1.28E+00 1.17E+00 1.23E+00 1.28E+00 1.31E+00 1.19E+00 

N/0 4.14E-01 4.20E-01 5.21E-01 4.70E-01 5.23E-01 5.67E-01 6.86E-01 

Ne/0 2.62E-01 2.65E-01 2.73E-01 2.68E-01 2.59E-01 2.71E-01 2.83E-01 



Table 6. Abundances determination for models with T* = 10 K. 



model Ab model Bb model Cb model Db model Eb model Fb model Gb 

T,[oin] 10388. 10588. 12968. 10974. 10597. 10472. 13067. 

T,[Nii] 10303. 10247. 10770. 10500. 10988. 11159. 11494. 

ne[sii] 3947. 3964. 4105. 3975. 2180. 5439. 6171. 

ne[Ariv] 5649. 5650. 5706. 5645. 8916. 8935. 9394. 



He+(A5876)/H+ 


1.07E- 


■01 


1.05E- 


■01 


8.67E^ 


■02 


1.05E-01 


1.05E^ 


■01 


1.08E-01 


8.58E^ 


■02 


He++(A4686)/H+ 


8.73E- 


■03 


I.IOE- 


■02 


2.69E^ 


■02 


1.13E-02 


9.09E^ 


■03 


8.96E-03 


2.82E^ 


■02 


He/H 


1.16E- 


■01 


1.16E- 


■01 


1.14E^ 


■01 


1.16E-01 


1.14E^ 


■01 


1.17E-01 


1.14E^ 


■01 


0+(A3727)/H+ 


5.70E- 


■05 


6.13E- 


■05 


7.34E^ 


■05 


6.02E-05 


3.09E- 


■05 


4.55E-05 


6.20E^ 


■05 


O++(A5007)/H+ 


3.79E- 


■04 


3.64E- 


■04 


2.86E^ 


■04 


3.65E-04 


4.08E^ 


■04 


3.51E-04 


2.61E^ 


■04 


0++(A1663)/H+ 


3.93E- 


■04 


3.79E- 


■04 


3.08E^ 


■04 


3.79E-04 


4.25E^ 


■04 


3.60E-04 


2.83E^ 


■04 


0++(A52^im)/H+ 


3.89E- 


■04 


3.76E- 


■04 


3.08E^ 


■04 


3.77E-04 


4.84E^ 


■04 


3.76E-04 


3.01E^ 


■04 


0++(A88^im)/H+ 


3.91E- 


■04 


3.77E- 


■04 


3.09E^ 


■04 


3.78E-04 


5.33E- 


■04 


4.33E-04 


3.52E^ 


■04 


0++(A4651)/H+ 


3.89E- 


■04 


3.76E- 


■04 


3.16E^ 


■04 


3.77E-04 


4.18E- 


■04 


3.67E-04 


2.94E^ 


■04 


O/H 


4.59E- 


■04 


4.55E- 


■04 


4.30E^ 


■04 


4.56E-04 


4.64E- 


■04 


4.18E-04 


3.91E^ 


■04 


C/O 


1.04E+00 


1.02E+00 


9.93E^ 


■01 


l.OOE+00 


1.03E+00 


1.06E+00 


9.97E^ 


■01 


N/O 


5.88E- 


■01 


5.80E- 


■01 


5.82E^ 


■01 


5.87E-01 


6.26E- 


■01 


8.12E-01 


7.55E^ 


■01 


Ne/O 


3.01E- 


■01 


3.04E- 


■01 


3.15E^ 


■01 


3.03E-01 


2.85E^ 


■01 


3.17E-01 


3.29E^ 


■01 



, the error in the derived O/H is much smaller. The reason is 
that photoelectric heating is strongest in the region populated 
by 0+++ ions, while the classical oxygen abundance determi- 
nation relies on lines emitted in the O^"*" and zones. For the 



model which is most affected by temperature inhomogeneities 
(the filamentary model with dust, i.e. model Gb), the Balmer 
jump temperature is 11 500 K. This is significantly lower than 
Ti[oiu] and than To(H+), as can be seen in Table I and For 
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such a nebula, the effect of dust heating would have dramatic 
consequences on the interpretation of ultraviolet lines from 
higher ionization species, such as and 

5.2. Interpretation of spatially resolved observations 

We have also computed the electron temperatures, densities and 
abundances for simulated observations on lines of sight located 
at different projected distances from the central star. The tem- 
perature profiles T^^q m] and Tj.[n ii] are shown in Fig. |^ for our 
dust-free homogeneous models (Aa and Ab) and for the homo- 
geneous models with small grains (Ca and Cb). We see that 
dusty models containing a population of small particles (our 
models Ca and Cb) do show a significant gradient in Tr[oiii] 
in contrast with dust-free models (Aa and Ab). For example 
Ti [Oiii] varies from 14200 K (in the center) to 10500 K (in the 
outskirts) in model Ca, and from 13500K (in the center) to 
8500 K (in the outskirts) in model Cb, while in dust-free mod- 
els the Trfoiii] variations are of the order of 500 K only. The 
variation of T,.[im uj across the face of the nebulae is very mild (a 
few hundred K only), even in the case of dusty models, which 
is not surprising since these lines arise in the outer parts of the 
nebulae. 



6. The temperature fluctuation parameter 

Numerous studies have used the temperature fluctuation con- 
cept introduced by Peimbert (1967) and derived P from obser- 
vational data (although as pointed out for example by Kingdon 
& Ferland 1995 the value derived from the observations does 
not strictly correspond to the definition of t^). Besides, as 
demonstrated by Stasinska (2001) on two-zone toy models, dif- 
ferent temperature distributions with formally the same value 
of and same mean temperature can give very different ratios 
of the various lines emitted by 

Because of this long tradition in the use of t^, it is worth 
commenting on the values of P found in our models (they cor- 
respond to the igjj. values of Kingdon & Ferland 1995). As seen 
in Tables 3 and 4, and as noted by previous authors, the values 
of t2(H+), i2(0+) and t^(0++) may actuafly be quite differ- 
ent. They are smaller than 10^^ for all our models except those 
containing small dust grains with a dust-to-hydrogen mass ra- 
tio 5 10^'^. Models with only one population of grains with the 
MRN size distribution between 0.01 and 0.25 /^m and a dust- 
to-hydrogen mass ratio of 10^^ (Models Ba and Bb) do give 
higher values of than the analogous dust-free counterparts, 
but they are still of the order of a few 10""^ or less. In this case, 
actually, higher values of result from a temperature gradient 
rather than a small scale temperature variation. The dust-free 
models with inhomogeneous density distributions return values 
of that are similar to constant density models. 

On the other hand, models containing small grains return 
values of largely above 10~^. In the constant density cases, 
this is due to the strong temperature gradients in the inner zone, 
which induce a i^(H+) of the order of 3 --4 10^^. For 0+, 
the value of P is smaller (but still around 10^^) because this 
zone is distant from the star and less affected by photoelectric 



heating. A similar remark holds for t^(0++) in the = 10^ K 
model. 

In the filamentary models with small grains (models Ga and 
Gb), The values of t^(H+) and <^(0++) are further increased 
with respect to the constant density models, because of the im- 
portant small-scale temperature variation seen in Fig. 7. This is 
the first time that photoionization models give such large val- 
ues of for nebulae with solar vicinity chemical composition 
and moderate effective temperatures. 

7. Concluding remarl<s 

In continuation of previous studie s on th e effects of dust grains 



in ionized plasmas (e .g. Spi tzer |1948| , Baldwin et al 



1991 



Weingartner & Draine |2001b[ ), we have investigated the role 
of dust in the thermal balance of planetary nebulae. This 
subject has been introduced by the works o f Bor kowski & 
Harrington ( [l 99 1| ) and Dopita & Sutherland ( ^00C| ), the first 
one studying the particular case of a hydrogen-poor dust heated 
planetary nebula, the second one exploring more generally the 
conditions under which dust heating plays a role in planetary 
nebulae. These last authors found that dust heating is impor- 
tant if a population of small dust grains is assumed. 

In the present study, we have analyzed in more detail the 
consequences of the presence of such grains in planetary nebu- 
lae, by considering several typical examples. In short, we have 
constructed photoionization models for two families of nebu- 
lae. One concerns rather dense nebulae (n(H) =10'' cm^-^) sur- 
rounding relatively cool central stars (T* =5 10'' K). The other 
concerns more diluted nebulae (n(H) = 5 lO'^cm^'') surround- 
ing hot central stars (T* = 10^ K). Such cases provide good ex- 
amples of different situations encountered in planetary nebu- 
lae. We have studied the temperature profiles of such models 
for different dust contents and different density distributions. 
As expected, the effect of dust heating is particularly important 
in the presence of a population of small grains. The electron 
temperature increases by a factor of about 2 in the central re- 
gions when adopting a population of small grains with a dust- 
to-hydrogen mass ratio of 5 10~^. Observational consequences 
of such a temperature increase are seen in optical spectra only if 
the high temperature zone contains a large proportion of 0++ 
ions which provide good temperature diagnostics. 

The presence of small grains in planetary nebulae can ex- 
plain a number of problems that have found no solution so 
far For example, tailored photoionization models of individ- 
ual planetary nebulae sometimes indicate an energy deficit in 
the thermal balance: the electron temperature derived from the 
observation of [O III] A4363/5007 is higher than can be ob- 
tained by photoionization models. This is the case in two of the 
five planetary nebulae with Wolf-Rayet central stars studied by 
Pena et al. ( |1998| ). 

Also, high signal-to-noise, spatially resolved observations 
of several planetary nebulae show that the temperatures derived 
from [O III] A4363/5007 smoothly decrease with projected dis- 
tance to the central star Indirect indications of such gradients 
are provided by ratios Trjomj/rriNii] larger than one in spa- 
tially resolved observations as well as in integrated spectra of 
many planetary nebulae. Such gradients are difficult to repro- 
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Fig. 8. Trfoiii] (thick lines) and T^jnh] (thin lines) at different projected distances from the central star for constant density 
dust-free models (Aa and Ab, left panels) and for constant density models with two populations of grains: small grains and large 
grains, each population having a dust-to-hydrogen mass ratio of 5 10^'^ (Models Ca and Cb, right panels). As in the other figures, 
upper panels correspond to T* = 5 10^ K and lower panels to T* = 10^ K. 



duce with photoionization models, even considering density 
gradients. They are naturally explained if one assumes that the 
nebulae contain small grains. 

Also, photoelectric heating by small dust grains naturally 
explains, at least qualitatively, why in planetary nebulae Balmer 
jump temperatures are smaller than the temperatures derived 
from [O III] A4363/5007. 

Detailed photoionization modelling of planetary nebulae 
often predicts too small a value of the [O i] A6300 intensity 
as compared to observations (Clegg et al. 1987, de Marco & 
Crowther 1999, Dudziak et al. 2000). One way to enhance this 
line in dust free plasmas is by assuming high density opti- 
cally thick condensations, which wiU strongly reduce the lo- 



cal ionization parameter and boost the [O l] A6300 lines with 
respect to smooth density models. However, [O l] A6300 can 
also be enhanced by additional heating in the transition zone 
between ionized and neutral gas. This heating can be provided 
by the photoelectric effect on small dust grains, whose impor- 
tance with respect to heating by photoionization increases as 
the Lyman continuum photons become exhausted. 

Moderate density inhomogeneities, such as are inferred 
from the distribution of the surface brightness of many plan- 
etary nebulae, produce only minimal spatial temperature vari- 
ations in dust-free nebulae. The presence of small dust grains 
completely modifies the small scale temperature structure of 
filamentary planetary nebulae, by spectacularly boosting the 
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temperature of the diluted component. This results in a situ- 
ation closer to pressure equilibrium between filament and am- 
bient gas than in absence of small grains. 

It is perhaps by measuring directly the temperature in and 
between filaments that one would have the best indirect indica- 
tion of the presence of small grains in planetary nebulae. 

Such observations would be really worthwhile, since our 
models show that small dust grains inside planetary nebulae are 
able to provide a natural solution to the long standing problem 
of temperature fluctuations first mentioned by Peimbert (1967). 

There is still a huge amount of work to be done on dust 
properties before photoelectric heating can be modelled to the 
same degree of reliability as other heating processes in ionized 
nebulae. Impressive progress in this sense has b een m ade re- 
centl y (We ingartner & Draine ZOOla, Li & Draine 2001 , Draine 
& Li 2001 , Weingartner & Draine 2001b ) but many uncertain- 
ties remain. 

We have shown that small dust grains are a promising ex- 
planation for a number of problems encountered in planetary 
nebulae. Still, several problems will not be solved. These are: 
i) the large inverse temperature gradients materialized by the 
small Ti.joiiij/IwNiii ratios observed in some planetary nebu- 
lae (Pena et al. 2001); ii) the very low Balmer jump temper- 
atures (around 3.5 10"^ K) observed in a few planetary nebu- 
lae (Liu et al. 2001a); iii) the origin of the large discrepancies 
between abundances derived from forbidden and collisionally 
excited lines of the same ions. Concerning the latter item how- 
ever, Liu et al. (2000) make the point that present computations 
of the effective recombination coefficients for line emission 
do not take into account dielectronic recombinations for states 
with principal quantum number above 10. These processes are 
expected to strongly enhance the emissivities of the recombina- 
tion lines at temperatures above 2 10^ K. Since small grains are 
able to boost the electron temperature to such high values in the 
most tenuous parts of filamentary or knotty planetary nebula, 
they might well also solve the recombination line conundrum, 
if the speculation about the enhancement of the recombination 
line emissivities at high temperatures proves to be true. 
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